Ceo: the first one-component gel? 
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Until now, gels have been formed of multicomponent soft matter systems, consisting of a solvent 
and one or more macromolecular or colloidal species. Here we show that, for sufhcient quench rates, 
the Girifalco model of Ceo can form gels which we identify by their slow dynamics and long-lived 
network structure. These gels are stable at room temperature, at least on the simulation timescale 
up to 100 ns. At moderate temperatures around 1000 K, below the bulk glass transition temperature. 
Ceo exhibits crystallisation and phase separation proceeds without the dynamical arrest associated 
with gelation, in contrast to many colloidal systems. 

PACS numbers; 



I. INTRODUCTION 



Gels form part of our everyday lives, yet are remark- 
ably poorly defined. At moderate timescales, gels lack a 
clear distinction from 'attractive glasses' at high densities 
[l[ , while at low densities gels share many characteristics 
of liquids. Recently however, a distinction has been made 
between equilibrium and non-equilibrium gels 0]. The 
former are associated with systems whose intermolecu- 
lar potential is not spherically symmetric, for example 
limited-valency models lead to equilibrium gels 0] and 
'empty liquids' Q. A second class of gels are formed 
through arrested metastable gas-liquid phase separation. 
In these 'spinodal gels' the metastable liquid phase is 
sufficiently dense to undergo dynamical arrest, suppress- 
ing demixing Q , an extreme case of viscoelastic phase 
separation [7|. It was recently shown that the forma- 
tion of clusters in this dense liquid phase prevents relax- 
ation to the underlying crystal Q . A third type of gel is 
formed through diffusion-limited cluster aggregation 0. 
Other types of gels form in polymer solutions (spinodal 
gels), precipitate from solution, for example silica, and 
are formed by clays such as laponite, the nature of which 
remains disputed [3. 

All these gels have one thing in common: they arc 
multi-component systems, consisting of a solvent and at 
least one macromolecular or colloidal component. Al- 
though in equilibrium, the degrees of freedom of the sol- 
vent and smaller macromolecular/coUoidal species can 
be formally integrated out, leading to a one-component 
treatment [ll|, it is important to note that considera- 
tion of an explicit solvent has recently been implicated 
in gelation, due to hydrodynamic interactions [l2|. 

Here we enquire whether gels can be formed from a 
single molecular species. It has emerged from the colloid- 
polymer mixture literature 0, that gels are associated 
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Figure 1: (color online) Interaction potentials. Solid black 
line, Girifalco potential used here to describe Ceo [IBl- Dashed 
red line, Asakura-Oosawa potential for colloid-polymer mix- 
tures with a size ratio of 0.18 [^. The Asakura-Oosawa pa- 
rameters are scaled to aid comparison with the Girifalco po- 
tential. 



with systems exhibiting a short-ranged attraction, com- 
pared to the particle/molecular diameter. Such systems 
also tend not to exhibit a stable liquid phase. Seeking 
a molecular system with such properties, one naturally 
turns to Ceo- While it is not yet clear whether Cgo has 
a stable liquid phase [l5l. the commonly used model 
potential of Girifalco [T] [l^, is understood to exhibit 
a liquid over a limited range in temperature (l7j . The 
short-ranged attractions in Cgo , relative to its molec- 
ular size, make it a most suitable candidate for one- 
component 'spinodal gelation'. We shall thus pose the 
question: can Cgo form a gel? 

Before proceeding, let us be clear about what we mean 
by gelation in this paper. Certainly, gels are percolat- 
ing networks. In the case of Cgo, dynamical arrest will 
be required to prevent phase separation on the measure- 
ment timescale. To create such a spinodal gel, we shall 
therefore quench Ceo hito the region of the phase dia- 
gram where gas-liquid phase separation may form a bi- 
continuous spinodal-like texture and undergo dynamical 
arrest. By dynamical arrest, we require that the perco- 
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Figure 2: (color online) Ground state clusters for CeoISJI- 
Here the clusters are defined according to our topological clus- 
ter classification. Ceo molecules are represented as spheres 
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Figure 3: (color online) Phase diagram for Ceo- Light dashed 
line is a guide to the eye denoting the liquid-gas binodal and 
solid line indicates freezing (from [TJl)- Percolation at low 
temperature is indicated by the solid black line. State points 
sampled with constant temperature simulations and corre- 
sponding structures are denoted as follows: (F) supercritical 
fluid, red circles; (F*) non-phase separated metastable fluid 
green circle; (X) phase-separated crystal, blue triangles; (C) 
isolated clusters, open square; gels, purple squares. 



lating structure persists on the simulation timescale. Our 
requirement that phase separation is suppressed on the 
molecular lengthscale distinguishes our approach from 
materials such as fumed silica, which is formed from 
a^ lomcratcs of particles of condensed amorphous silica 

As in a previous study on a colloidal system Q, 
we shall characterise the structure in terms of clusters 
formed by small groups of Cgo molecules in isolation. The 



structures formed are then analysed in terms of topolo- 
gies based on the minimum energy clusters using an al- 
gorithm we have recently developed for identifying clus- 
ters in condensed phases, the topological cluster clas- 
sification (TCC) [13 • Clusters of Cgo molecules have 
received considerable attention [2^ and for small clus- 
ters the ground states are known both experimentally 
[20l - [23j and computationally 24 1. Moreover, numerical 
studies have shown that for larger clusters Ceo exhibits 
kinetic trapping, leading to a discrepancy between pre- 
dicted ground state structures [l^l and those observed 
experimentally below 500 K, which required annealing 
to access the ground state [22|, [2^ . 

We shall use molecular dynamics simulation to capture 
the behaviour of Cgo on a timescale and system size suf- 
ficient that we may discuss gelation. For this purpose, 
we use the intermolecular interaction potential proposed 
by Girifalco [l^, which reads 



u{s) = ai 
+a2 



s(s- 1)3 s(.s + l)3 s4 
1 1 2 

s(s- 1)9 ^ s(s + l)9 " 7^ 



(1) 



r/2a and 2a 



7.1 Angstroms, 

'18 



ai = 74.94 X 



where s 

10-15 erg, a2 = 135.95 x lO"^* erg. The Girifalco po- 
tential is plotted in [TJ This Girifalco model is known 
to undergo dynamical arrest at high temperature and 
density j25l - [27| . The glass line intersects the gas-liquid 
spinodal at around To ~ 1100 K, while crystallisation is 
expected at higher temperature [23. While locally crys- 
talline structures still meet our criterion for gelation (ie 
a long-lived network), the absence of vitrification of the 
dense phase suggests that demixing may proceed. To 
prevent demixing, we need to quench below Tq . Noting 
that free surfaces enhance local diffusion [1, [2^ [2^ , we 
expect that it will be necessary to quench well below Tq. 
This is consistent with experimental work where anneal- 
ing to 610 K enabled small clusters to equilibrate [23j . 
Figure 1 also depicts the Asakura-Oosawa (AO) poten- 
tial [13, HH often used to described colloid-polymer mix- 
tures. The AO potential is plotted here to correspond to 
a polymer-colloid size ratio of 0.18 Q which leads to a 
relative interaction range rather shorter than that of the 
Girifalco potential. 

Our procedure is therefore as follows: we simulate Gir- 
ifalco Ceo at constant temperature, following an 'instan- 
taneous' quench to enquire at what, if any, temperature. 
Ceo forms a gel stable on the simulation timescale. We 
then perform quenches from high temperature to deter- 
mine what quench rate is required to form gels. 
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Figure 4: (color online) Typical (non-equilibrium) structures, (a) non- phase separated fluid ( p = 0.43 nm"^, T = 1400 K), (b) 
phase separated crystal (p = 0.430 nm"^, T = 1000 K), (c) gel ( p = 0.215 nm"^, T = 300 K); (d) isolated clusters (p = 0.108 
nm"'^, T = 300 K). All data shown here are from simulations at constant temperature. Colors denote the differing clusters or 
crystal structures to which the molecules belong, as shown on the right. Free denotes molecules not belonging to any clusters. 
Molecules are drawn around half actual size and represented as spheres. 



II. SIMULATION DETAILS 

A. Molecular dynamics simulation 

Wc carry out simulations of two types. Constant tem- 
perature simulations, where the system is run up to 
10^ timesteps (99.6 ns) and sampled for a further 10^ 
timesteps, and temperature quenches where temperature 
is decreased from 2328 K to 300 K at rates between 
2.04 X 10" and 2.04 x 10" Ks'^ State points sam- 
pled for constant temperature simulations are shown in 
|3l Throughout we use molecular dynamics simulation 
with a Gaussian thermostat [H] with a timestep of 9.96 
fs. We typically used 4000 or 2048 Cgo molecules. Runs 
with 4000 particles showed no qualitative differences. Oc- 
casionally, long times required smaller simulations of 864 
particles. Data are shown after 9.96 ns of 'equilibration' 
or 'aging', unless otherwise stated. All simulations are 
performed at constant volume. We focus on data taken 
along the critical isochore p ~ 0.43 nm~^ unless other- 
wise stated. Results at lower density p ~ 0.215 nm^^ are 
qualitatively similar. 

B. The topological cluster classification 

To analyse the structure, we identify the bond network. 
Here we define two molecules as bonded if they approach 
within 1.23 nm, which is the approximate location of the 
first minimum of the pair correlation function g{r) upon 
condensation. Having identified the bond network, we 
use the Topological Cluster Classification (TCC) to de- 
termine the nature of the cluster [l^ . This analysis iden- 
tifies all the shortest path three, four and five membered 
rings in the bond network. We use the TCC to iden- 
tify structures of between 5 and 13 molecules which are 
topologically equivalent to ground states for the Girifalco 
potential ([T]) [2J]. The clusters we identify are depicted 
in In addition we identify the FCC and HCP thir- 
teen particle structures in terms of a central particle and 



its twelve nearest neighbours. For spherically symmetric 
potentials differing cluster structures are found as a func- 
tion of interaction range [s^. For 11-membered clusters, 
our topological approach does not distinguish between 
Cs (IID in [13) and C2v (HC in [H). Likewise for 12- 
membcrcd clusters Cg (12C in [s^) and C^y (12B in psl ) 
are identical for our purposes. For more details see [l9[. 
If a molecule is found to be part of more than one cluster, 
we count it as the larger cluster. Moreover, if a molecule 
is found to be part of both a HCP and FCC local envi- 
ronment (for example in random close packed stacking), 
we count it as FCC. 



III. RESULTS 

A. Constant temperature simulations 

We begin our presentation of the results by discussing 
our interpretation of the different structures we find. The 
equilibrium phase diagram is sketched in [3] [l7| and typ- 
ical structures are shown in 21 However, the system 
is not always able to fully equilibrate on the simula- 
tion timescale, a point to which we return below. For 
the densities we consider, at temperatures less than the 
triple point (around 1880 K [13), equilibrium for the Gir- 
ifalco model of Ceo is gas-crystal coexistence. However, 
for weak degrees of undercooling, crystal nucleation may 
not occur on the system sizes and timescales we con- 
sider, while for deep quenches, dynamical arrest limits 
access to equilibrium. In our simulations, these two non- 
equilibrium states are manifested as a metastable fluid, 
which can (but does not always) undergo phase separa- 
tion prior to crystallisation, and a gel which forms a per- 
colating network respectively. We identify these states by 
visual inspection of coordinate data and also require that 
more that 50% of molecules are in a locally HCP or FCC 
environment according to the TCC for any state point to 
be considered to be crystalline. We further require that 
gel states exhibit slow dynamics, in other words, that 
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Figure 5: (color online) Mean squared displacement (a) and 
pair correlation functions (b) for constant temperature simu- 
lations. Here p = 0.43 

nm In (b) the data are offset for clarity. 



they do not phase separate or significantly coarsen on 
the simulation timescale. 

More rigorous approaches, notably testing for percola- 
tion, have the disadvantage that for the system sizes we 
are able to tackle here, for p ^ 0.2 nm~^, all state points 
were found to percolate according to our bonding crite- 
rion, despite the different structures evident in|4l At low 
density, isolated clusters were found [4|d)]. In our simu- 
lations, these clearly do not form a percolating network. 
Thus we sketch a percolation line in [31 

We have argued that in the case of Cgo, gelation should 
be associated with slow dynamics. Plotting the mean 
squared displacement (MSD) at different temperatures 
inlSfa), we see a sharp fall in MSD between T=1400 and 
1000 K. Thus we expect that, at 1000 K and below, we 
may find gelation. Inspection of the coordinate data in 
mjb) shows that at 1000 K, Ceo crystallises on our simula- 
tion timescales. This is consistent with nucleation studies 
carried out at higher density (s^ . However at some of 
the lower temperatures we consider, it is clear Cgo is suf- 
ficiently arrested that little restructuring is seen on the 
simulation timescale. Dynamically, therefore, it is possi- 
ble that Ceo can satisfy our requirements for classifica- 
tion as a gel. We note that these measures are performed 
on a system which is not in equilibrium for some state 
points. However, our intention here is to qualitatively 
demonstrate a substantial slowing of the dynamics. 

Turning to local structural measures, namely the pair 
correlation function g{r) we see three regimes injSjb). At 
high temperatures, a fluid g{r) is found. Note that even 
at temperatures below the triple point |17| . we see little 
change in g{r) on these simulation timescales. However, 
at lower temperatures (1000 K), the system crystallises. 
At lower temperatures still (300 K), although there is 
some local structure, overall the ordering is reduced. For 
temperatures less than around 1000 K g{r) takes values 
substantially greater than unity at short range. This in- 
dicates density fluctuations on the lengthscale of a few 
particle diameters. While the limited size of our simula- 
tions preclude a definitive statement, we note that such 
local density fluctuations are indicative of networks, such 
a that found in (UJc) . 
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Figure 6: (color online) The population of molecules in ground 
state clusters of Ceo as a function of inverse temperature. 
Here p = 0.43 nm~'^. Shaded regions are a guide to the eye. 
Note the semi-log scale. 
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Figure 7: (color online) (a) Static structure factor S(q) for 
gels at different ageing times t^. Black line tn, ~ 1 ns, pink 
line, = 10 ns, red line, = 100 ns. Inset shows first peak. 
p = 0.21 

nm^'^ T — 300 K. (b) Mean squared displacement for surface 
(pink), and interior (purple) molecules, p = 0.43 nm~^ 
T = 300 K. 



Using the topological cluster classification, the same 
three regimes are identified, as shown in|ni At high tem- 
peratures (T>1400 K), the population of all clusters rises 
sharply upon reducing the temperature. Around lOOOK, 
the structure is dominated by molecules in a locally crys- 
talline environment, and the population of other clusters 
drops. At lower temperatures still, the population of lo- 
cally FCC and HCP molecules drops, and the system is 
once again dominated by amorphous clusters. Of these, 
at 300 K, 11-mcmbcred Cg clusters arc the most numer- 
ous, accounting for some 20% of the total population of 
Ceo molecules. These clusters are five-fold symmetric, 
and it is tempting to infer that this five-fold symmetry 
might impede crystallisation (sBj . The five- fold symmet- 
ric structure popularised by Frank (35| . the IS-membered 
icosahedron is less popular than either HCP or FCC in 
this case, despite being the minimum energy cluster for 
13 molecules [2J]. We note that the icosahedron has been 
found to be relatively resistant to kinetic trapping [s^ 



and is also observed experimentally [21 



We now consider the gel stability. Spinodal gels such 



5 




((ns) ?-(nm) 

Figure 8: (color online) Mean squared displacement (a) and 
pair correlation functions (b) for quenches. Here p — 0.43 
nm~^. Quench rate is denoted in Ks^^. 

as those we consider here are metastable to phase sepa- 
ration. Although that process is suppressed by slow dy- 
namics induced by quenching, it is reasonable to pose the 
question, to what extent is phase separation effectively 
suppressed? To this end, wc plot the static structure fac- 
tor S{q) in Fig. [Tl^a). Indeed a small degree of structural 
evolution is found, as shown by the slight shift in the 
first peak of S{q). While this indicates some degree of 
coarsening, it is very clear that the gel lifetime exceeds 
the simulation timescale, noting that here the simulations 
are run for up to 100 ns. 

Nevertheless, there is some degree of coarsening. We 
investigate possible coarsening mechanisms by consider- 
ing the mean squared displacement of particles on the 
surface and in the interior of the 'arms' of the gel. We 
define molecules on the surface as those with fewer than 
eight neighbours. The analysis in Fig. [TJb) shows that 
indeed, molecules on the surface have a marginally higher 
mobility than those in the interior. This is qualitatively 
consistent with previous studies of gelation in colloidal 
systems [i,[2^. 



B. Temperature quenches 

The picture that emerges from our analysis of the con- 
stant temperature runs is that at T<600 K, Ceo can form 
gels, according to our criteria. However such 'instanta- 
neous quenches' are unphysical for molecular systems (al- 
though essentially accessible to soft matter systems such 
as colloids). We therefore enquire as to what quench rate 
is required to form a gel. We choose to quench from a 
supercritical fluid at 2328 K to 300 K. Inspection of the 
phase diagram suggests that, for sufficiently slow quench 
rates, we might expect phase separation to occur in tem- 
perature regime 1940 K > T>600 K, and that crystalli- 
sation may also occur. 

We thus consider a range of quench rates, and identify 
the final state as before. Data are presented for the 9.96 
ns immediately following the quench. Mean squared dis- 
placement data [HJa)] shows that the quenches we have 
performed result in a state with relatively slow dynamics. 
Mobility is strongly suppressed for slower quench rates. 



The pair correlation function data in ID^b) show a sim- 
ilar structural transitions to the constant temperature 
data, although here the transition occurs as a function 
of quench rate. For quench rates at or above 2.04 x 10^^ 
Ks^^, we find that the g{r) remains amorphous. A slower 
quench rate of 2.04 x 10^^ Ks~^rcsults in a crystalline 
structure, which undergoes phase separation, similar to 
Hb). 

This observation is confirmed with the TCC analy- 
sis, where we consider the cluster populations in the 
quenched systems. In[5Ja) wc show average cluster popu- 
lations for runs of 9.96 ns following the quench, and find a 
clear shift to a crystal-dominated structure for the slow- 
est quench rate. In fact the trend towards a higher de- 
gree of crystalinity is apparent between the higher quench 
rates of 2.04 x lO^^ ^nd 2.04 x lO^^ ^^-i Turning to the 
mean squared displacement data in[5Ka) we note that for 
rapid quenching in particular, some movement persists. 
In other words, the gel undergoes aging. Restructuring 
associated with this aging is apparent in the larger error 
bars for the 2.04 x lO^^ Ks-^quench rate inlKa). InMh) 
we consider the local structural consequences of this ag- 
ing. The data shows a gradual trend of increasing HCP 
and FCC local structures and a decrease of the most pop- 
ular cluster, m ~ 11 C2v This increase in crystalinity 
is consistent with previous studies on aging colloidal gels 



IV. DISCUSSION AND CONCLUSIONS 

We begin the discussion by considering how practi- 
cally realisable the Ceo gels we have predicted might be. 
We have shown that gels form and exist on a simulation 
timescale of 10 ns. This leads to two main questions: 
firstly, can the gels form in the first place and secondly, 
how long can they last? Evidence from experimental 
work on clusters of Ceo molecules suggests that in fact, 
long-lived metastable states are possible. Below around 
500 K, clusters show little preference for ground states 
[20I . [23 . . Wc thus argue that Ceo can be prevented 
from equilibrating on an experimental timescale at room 
temperature, although we note that some aging can be 
seen in our simulations. The first question is potentially 
more challenging. The quench rates we have used are 
not currently experimentally feasible (maximum quench 
rates are around 10^ Ks~^by vapor deposition). Our pre- 
dictions arc that Cgo will phase separate to crystal-gas 
coexistence on such a timescale. However, shear could be 
used to prevent crystallisation during quenching, allow- 
ing milder quench rates. Other possibilities include com- 
pressing low-temperature amorphous clusters. We note 
in m^d) that the isolated clusters formed at T= 300 K 
are far from spherical and might reasonably be expected 
to form a gel upon increasing the density. Following [s^ 
we also note that the Girifalco potential is not the only 
model for Ceo- The Pachcco-Prates-Ramalho model [S^ 
is in fact more sticky (shorter ranged) than that of Giri- 
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Figure 9: (color online) The population of molecules in ground 
state clusters of Ceo for different quench rates (a). Quench 
rates are denoted in Ks^^. (b) Structural changes upon aging 
in a gel formed after a quench rate of 2.04 x lO^^Ks^^. Here 
p = 0.43 nm~^. 



falco, which would be expected to promote gelation. 

Another approach to producing one-component gels 
might be to consider larger fuUerenes whose dynamics 
will be much slower. Such a candidate is C54owhich may 
also exhibit a molecular shape with high symmetry and 
stability [39| , although its high binding energy would ne- 
cessitate quenching from very high temperatures at which 
thermal stability of the fullercne molecule may become a 
problem. The same might even hold for Cgo whose ther- 
mal stability at high temperatures (such as criticality) is 



disputed I40|,|4l|, although at temperatures below 1700 
K Cgo is likely stable, so compressing disordered clusters 
remains a potential experimental route to forming gels of 
Ceo- 

The behaviour exhibited by Cgo in the constant tem- 
perature simulations may be compared with that ob- 
served in colloidal suspensions, where interactions may 
be tuned [l^, and which typically undergo 'instanta- 
neous' quenching, due to the slow relaxation time of the 
mesoscopic colloidal particles. The main difference be- 
tween the results presented here and those from a pre- 
vious study on a colloidal gel [1] are that Ceo exhibits 
a regime of crystallisation, which is almost entirely ab- 
sent from the colloidal system, where the crystal-like lo- 
cal structures accounted for no more than one particle in 
1000. There are at least three possible reasons for this 
difference. The first is that the relative interaction range 
is substantially longer in the case of Ceo and this may 
facilitate crystallisation (see [T|). Secondly, colloids are 
polydisperse, which might suppress, although should not 
prevent, crystallisation [Z^, E^- Finally, we have used 
Newtonian dynamics for Ceo while colloids are diffusive, 
although we are unaware of any studies considering the 
role of dynamics in crystallisation. In colloid simulation 
work with Brownian dynamics with parameters com- 
parable to experiments Q there is some evidence that 
crystallisation can occur on the low temperature side 
of the metastable gas-liquid binodal in a monodisperse 
system. This suggests that the role of polydispcrsity in 
the colloid experiments may be important in suppressing 
crystallisation. 

In summary, we have presented numerical evidence 
that Ceo, under the right conditions can form a gel, which 
forms through arrested spinodal decomposition. Poten- 
tially, therefore, one-component gels might be realised. 
Whether these conditions are experimentally accessible 
remains an open question. 
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